i 


AD-A243  288 


•  •  » 


:CH  Laboratory  of  Electronics 

TS  INSTITUTE  OF  TECHNOLOGY 
CAMBRIDGE,  MASSACHUSETTS  02139-4307 


November  15,  1991 


Dr.  Charles  A.  Luther,  Code  1121RS 
Scientific  Officer 

Office  of  the  Chief  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 


Re:  Contract  NOOOm-89- J- 1 1  07 


Type  of  Material; 

Title: 

Submitted  By: 

Period  Covered: 

No.  of  copies: 

Distribution : 


Annual  Technical  Report 

Active  and  Passive  Remote  Sensing  of  Ice 

J.A.  Kong 

October  1,  1990  -  September  30,  1991 

Three  plus  form  298 

Navy  Distribution  List  (9) 


Mary  S.  Greene 

RLE  Publications  Office 

Room  36-*4l2 


cc ;  J.A.  Kong  ( 1  ) 

D.  Duffy 

A.F.  Favaloro,  E19-702 
File  (1) 


OSP  71367 


DISTRIBUnON  STA'TBMEjrr  A 

Approved  lor  public  rdeoie; 
Dletributtoa  UnBmited 


91-16073 

illlllBH 


n‘>o 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OM$  No.  0704-018$ 


HMc  ftoentng  tertftA  for  mi%  colitaion  of  mformatiert  •»  «tcim«tto  to  «wor»9»  1  tiour  ptr  foipome.  tfidoding  the  time  for  rovtowiftq  imtrucoom,  Morchino  eiistmg  Oau  io^'A 
- - - - - - - - - -  - ~ ~  nd f»>i»wmg tho coi*oct>oo of  mformtioo.  trrrf  rnri-t*n  “T|T‘if*-^ tliti  tnnft*  tittmitf  nr  ini  rtrhtr  iifurr  ~~  ~n  i 


OoM  Highway.  Sort*  120*.  AfHogfon.  v*  222024302.  aod  10  tht  Offk*  of  MoAogofi^t  miO  Oudgot.  faptrwork  Roduction  Protect  P70*<01M).  Woitim9toA.  OC  20903. 


1.  AGENCY  USE  ONLY  (Ltst^  bUnk)  12.  REPORT  DATE  13.  REPORT  TYPE  AND  OATES  COVERED 

11-15-91  Annual  10/1/90  - 


4.  TITLf  ANO  SUBTI7U  5.  FUNDING  NUMBERS 

Active  and  Passive  Remote  Sensing  of  Ice 


C.  AUTHOR(S) 


Prof.  J.A.  Kong 


N00014-89-J-1 107 


425F023-03 


7.  KRFORMMIG  ORGANIZATION  NAME(S)  ANO  AOORESS(ES} 

Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS<ES) 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  In  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


IZa.  DISTRIBUTION /AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  200  words) 


Work  by  Prof.  J.A.  Kong  and  his  collaborators  is  summarized  here 


IS.  NUMBER  OF  PAGES 


IS.  PRICE  CODE 


17.  SECURITY  CLASSIFIUTION  111.  SECURITY  CUSStFICATION 
OF  REPORT  I  OF  THIS  PAGI 


IB.  SECURITY  CLASSIFICATIONI 20.  UMITATION  OF  ABSTRACT 
OF  ABSTRAa  I 


unclassified 


NSN  7540.01 -280-SS00 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rav.  2>B9i 
.raicnMa  »r  **«•  M  ZW-lt 
m-iu 


ANNUAL  REPOBT 


Title:  ACTIVE  AND  PASSIVE  REMOTE  SENSING  OF  ICE 
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Office  of  Naval  Research 

Contract  number:  N00014-89-J-1107 


ACTIVE  AND  PASSIVE  REMOTE  SENSING  OF  ICE 


ANNUAL  REPORT 


This  annual  report  covers  research  under  the  sponsorship  of  the  ONR  contract 
N00014-89-J-1107  from  October  1,  1990  to  September  30,  1991.  We  have  published  15 
journal  and  conference  papers  and  3  student  theses  during  this  period. 


1.  Multi-Layer  Random  Medium  Model  for  Sea  Ice 

Fully  polarimetric  scattering  of  electromagnetic  waves  from  snow  and  ice  is  stud¬ 
ied  with  a  multi-layered  random  medium  model  and  applied  to  interpret  experimental 
data  obtained  under  laboratory  controlled  conditions  such  as  CRRELEX.  The  snow  layer 
is  modeled  as  an  isotropic  random  medium.  The  sea  ice  is  described  as  an  anisotropic 
random  mediiun  due  to  the  nonspherical  shape  of  brine  inclusions.  The  underlying  sea 
water  is  considered  as  a  homogeneous  half-space.  The  random  media  in  both  layers  are 
characterized  by  three-dimensional  correlation  functions  with  variances  and  correlation 
lengths  corresponding  to  the  fluctuation  strengths  and  the  physical  geometries  of  the  inho¬ 
mogeneities,  respectively.  The  strong  fluctuation  theory  is  used  to  calculate  the  effective 
permittivities  of  the  random  media.  The  distorted  Bom  approximation  is  then  employed 
to  obtain  the  covariance  matrix  which  represents  the  fully  polarimetric  scattering  prop¬ 
erties  of  the  snow-ice  media.  It  has  been  shown  that  the  polarimetric  covariance  matrix 
contains  more  information  than  the  conventional  scattering  coefficients  on  the  remotely 
sensed  media. 


In  saline  ice  sheets  under  quiescent  condition,  the  background  ice  grows  in  colum¬ 
nar  form  and  saline  water  is  trapped  between  ice  platelets  in  the  form  of  brine  inclusions 
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which  are  usually  ellipsoidal.  The  ice  tends  to  grow  vertically  downward  rendering  the 
ellipsoidal  inclusions  aligned  preferably  in  the  vertical  direction  and  the  crystallographic 
C  axes  parallel  to  the  horizontal  plane.  In  this  case,  the  C  axes  are,  however,  random  in 
azinrathal  direction.  The  strong  fluctuation  theory  is  extended  to  account  for  vertically 
aligned  ellipsoidal  brine  inclusions  with  C  axes  randomly  oriented  in  the  horizontal  direc¬ 
tion.  The  brine  inclusions  are  described  by  three-dimensional  local  correlation  functions. 
The  configuration  average  over  the  azimuthal  orientation  angles  is  carried  out  in  the  pro¬ 
cess  of  deriving  the  global  correlation  tensor.  The  distorted  Bom  approximation  is  applied 
to  obtain  the  covariance  matrix  for  the  multi-layered  snow-ice  configuration.  The  theo¬ 
retical  results  show  non-zero  cross-polarized  returns  under  the  first-order  distorted  Bom 
approximation.  We  have  also  compared  the  results  with  experimental  data  obtained  by 
the  US  Army  Cold  Regions  Research  and  En^neering  Laboratory  (CRREL). 

In  sea  ice,  the  scatterers  can  have  various  shapes,  sizes,  and  permittivities.  We 
have  also  investigated  the  modeling  of  radar  backscatter  from  random  media  with  multiple 
scatterer  species.  We  consider  each  type  of  scatterers  as  a  species  which  can  take  on  a  shape, 
size,  and  complex  permittivity  difierent  from  other  species.  The  multiple  species  in  the 
random  medium  are  considered  as  randomly  oriented  ellipsoids  and  described  by  multiple 
three-dimensional  ellipsoidal  local  correlation  functions.  The  variances  and  correlation 
lengths  of  the  correlation  functions  characterize  the  fluctuation  strengths  and  the  pliysical 
geometries  of  each  species  of  scatterers.  The  effective  permittivity  of  the  random  medium 
is  derived  under  the  strong  fluctuation  theory  and  the  polarimetric  scattering  coefficients 
are  calculated  for  the  layer  configuration  with  the  distorted  Bom  approximation.  Due 
to  the  non-spherical  shape  and  the  random  orientation  of  the  scatterers,  the  correlation 
coefficient  between  the  HH  and  W  returns  has  a  magnitude  different  from  unity  and  a 
small  phase  angle.  The  scattering  coeffidents  are  also  used  to  calculate  the  Mueller  matrix 
for  synthesis  of  polarization  signatures.  The  co-polarized  signature  of  the  random  medium 
has  a  rather  straight  distortion  track  and  a  recognizable  pedestal. 


2.  Polarimetric  Radar  Calibration 


Accurate  calibration  of  polarimetric  radar  systems  is  essential  for  the  polarimetric 
remote  sensing  of  earth  terrain.  Polarimetric  calibration  algorithms  using  in-scene  re¬ 
flectors  and  distributed  targets  are  studied.  The  transmitting  and  receiving  ports  of  the 
polarimetric  radar  are  modeled  by  two  unknown  polarization  transfer  matrices.  These 
unknown  matrices  are  determined  using  the  the  measured  scattering  matrices  from  the 
calibration  targets.  For  the  case  of  polarimetric  calibration  using  three  in-scene  reflec¬ 
tors,  a  Polarization-basis  IVansformation  (PT)  technique  is  introduced  to  And  out  a  new 
transmitting  and  receiving  polarization  basis  under  which  the  scattering  matrices  of  the 
calibration  targets  will  fall  into  one  of  six  simpler  sets.  The  calibration  solution  can  then 
be  solved  easily  in  the  new  polarization  basis  and  converted  to  obtain  the  solution  in  the 
original  polarization  basis.  The  uniqueness  of  polarimetric  calibration  using  three  targets 
is  addressed  for  all  possible  target  combinations.  The  PT  technique  can  also  be  applied 
to  the  polarimetric  calibration  using  a  combination  of  arbitrary  in-scene  reflectors  and 
distributed  targets.  The  eflfect  of  misalignment  of  calibration  targets  and  the  sensitivity 
of  polarimetric  calibration  algorithms  to  the  noise  are  illustrated  by  investigating  several 
sets  of  calibration  targets. 

S.  Classiflcation  Studies  of  Polarimetric  SAR  Images 

We  have  studied  the  SAR  image  classification  by  using  the  neural  network  meth¬ 
ods.  Supervised  methods,  including  both  conventional  Maximum  Likelihood  (ML)  and 
more  recent  multi-layer  perceptron  neural  network  classifiers  have  yielded  higher  accuracy 
than  unsupervised  techniques,  but  suffer  from  the  need  for  human  interaction  to  prede¬ 
termine  classes  and  training  re^ons.  In  contrast,  unsupervised  methods  determine  classes 
automatically,  but  generally  show  limited  ability  to  accurately  divide  terrain  into  natural 
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dasses.  We  introduced  a  new  terrain  dassification  technique  to  discriminate  sea  ice  sig¬ 
natures  in  polarimetric  SAR  images  by  utilizing  imsuperrised  neural  networks  to  provide 
automatic  dassification,  but  employing  an  iterative  algorithm  which  overcomes  the  poor 
accuracy  of  other  unsupervised  techniques. 

Several  types  of  unsupervised  neural  networks  are  first  applied  to  the  dassification 
of  SAR  images,  and  the  results  are  compared  with  those  of  more  conventional  unsupervised 
methods.  Neural  network  approaches  indude  Adaptive  Resonance  theory  (ART),  Learning 
Vector  Quantization  (LVQ),  and  Kohonen’s  self-organizing  feature  map.  Conventional 
dassifiers  utilized  are  the  migrating  means  dustering  algorithm  and  the  K-means  clustering 
method.  Preprocessing  is  performed  with  the  SAR  images  to  reduce  speckle  noise  and 
stabilize  the  training  process  for  the  neural  networks.  Results  after  preprocessing  show  that 
LVQ  and  Kohonen’s  self-organizing  feature  map  outperform  the  conventional  unsupervised 
dassifiers,  but  are  still  inferior  to  supervised  methods. 

To  overcome  this  poor  accuracy,  an  iterative  algorithm  is  constructed  where  the 
SAR  image  is  redassified  using  a  Maximum  Likelihood  (ML)  classifier.  Training  of  the  ML 
dassifier  is  performed  using  a  training  data  set  first  classified  by  the  above  unsupervised 
method,  thus,  requiring  no  human  intervention,  and  preserving  the  unsupervised  nature 
of  the  overall  classification  scheme.  The  process  is  then  repeated  iteratively,  training  a 
second  ML  dassifier  using  data  dassified  by  the  first.  It  is  shown  that  this  algorithm 
converges  rapidly,  and  significantly  improves  dassification  accuracy.  Performance  after 
convergence  is  seen  to  be  comparable  to  that  obtained  with  a  supervised  ML  dassifier, 
while  maintaining  the  advantages  of  an  unsupervised  technique. 

The  unsupervised  and  iterative  techniques  devdoped  have  been  applied  to  the 
polarimetric  SAR  images  of  Beaufort  sea  ice  acquired  by  the  C-,  L-,  and  P-band  SAR 
instruments  of  Jet  Propulsion  Laboratory.  The  results  obtained  with  the  new  algorithms 


are  compared  with  the  results  obtained  with  other  techniques  by  classifying  terrain  features 
in  polarimetric  SAR  images. 

4.  K-Distribution  and  Polarimetric  Terrain  Radar  Clutter 

A  multivariate  K-distribution  has  been  developed  to  model  the  statistics  of  fully 
polarimetric  radar  data  from  earth  terrain  with  polarizations  HE,  HV,  VH,  and  VY.  In 
this  approach,  correlated  polarizations  of  radar  signals,  as  characterized  by  a  covariance 
matrix,  are  treated  as  the  sum  of  N  n-dimensional  random  vectors;  N  obeys  the  negative 
binomial  distribution  with  a  parameter  a  and  mean  N.  Subsequently,  an  n-dimensional 
K-distribution,  with  either  zero  or  nonzero  mean,  is  developed  in  the  limit  of  infinite  N  or 
illuminated  area.  The  probability  density  function  (PDF)  of  the  K-distributed  vector  nor¬ 
malized  by  its  Euclidean  norm  is  independent  of  the  parameter  a  and  is  the  same  as  that 
derived  from  a  zero-mean  Gaussian-distributed  random  vector.  The  above  model  is  well 
supported  by  experimental  data  provided  by  MIT  Lincoln  Laboratory  and  the  Jet  Propul¬ 
sion  Laboratory  in  the  form  of  polarimetric  measurements.  The  results  are  illustrated  by 
comparing  the  higher-order  normalized  intensity  moments  and  cumulative  density  func¬ 
tions  (CDF)  of  the  experimental  data  with  theoretical  results  of  the  K-distribution. 

5.  Composite  Volome  and  Surface  Scattering  Model 

Among  the  various  theoretical  models  applied  to  study  the  electromagnetic  wave 
scatterings  from  geophyncal  terrain,  such  as  snow  and  ice,  the  radiative  transfer  theory 
has  drawn  intensive  attention  in  the  microwave  remote  sensing  society  during  the  past 
years.  In  most  of  the  scattering  models,  the  volume  scattering  and  the  surface  scattering 
effects  have  been  investigated  separately.  Recently,  there  has  been  a  growing  interest  in 
the  construction  of  composite  models  which  can  take  into  account  both  types  of  scattering. 


We  derived  the  first  order  iterative  solution  to  the  vector  radiative  transfer  equations  for  a 
two-layer  medium  with  a  diffuse  top  boundary  and  an  irregular  bottom  boundary  of  Gaus¬ 
sian  roughness.  The  Kirchhoff  approximation  and  the  geometrical  optics  approach  with 
shadowing  correction  are  used  in  formulating  the  boundary  conditions.  To  demonstrate 
the  utilities  of  the  theory,  randomly  oriented  spheroidal  discrete  scatterer  model  is  used  to 
calculate  the  backscattering  coefficients  from  soybean  field  in  different  growing  stages  and 
compared  to  the  experimental  measurements.  Good  agreement  has  been  achieved  for  both 
the  co-polarized  and  the  cross-polarized  data.  It  is  observed  that  the  presence  of  the  rough 
surface  can  significantly  enhance  the  backscattering  at  small  incident  angles  and  increase 
the  cross-polarized  returns.  The  polarization  signatures  calculated  based  on  the  Mueller 
matrix  show  a  straight  distortion  track  and  an  observable  pedestal.  Numerical  comparison 
to  the  backscattering  coefficients  calculated  by  using  plainar  bottom  boundary  conditions 
with  or  without  the  incoherent  siddition  of  the  rough  stirface  effects  are  ailso  made. 

6.  Polarimetric  Passive  Remote  Sensing 

The  concept  of  polarimetry  in  active  remote  sensing  is  extended  to  passive  remote 
sensing.  The  potential  use  of  the  third  and  fourth  Stokes  parameters  U  and  V,  which 
play  an  important  role  in  polarimetric  active  remote  sensing,  is  demonstrated  for  passive 
remote  sensing.  It  is  shown  that,  by  the  use  of  the  reciprocity  principle,  the  polarimetric 
parameters  of  passive  remote  sensing  can  be  obtained  through  the  solution  of  the  associated 
direct  scattering  problem.  In  particular,  the  full  polarimetric  information,  including  the 
corresponding  brightness  temperatures  of  U  and  V,  can  be  obtained  from  the  solution  of 
the  direct  scattering  problem  for  four  different  polarizations  of  the  incident  wave.  These 
ideas  are  applied  to  study  polarimetric  passive  remote  sensing  of  periodic  surfaces.  The 
solution  of  the  direct  scattering  problem  is  obtained  by  an  integral  equation  formulation. 
Incidence  on  a  penetrable,  lossy,  medium  is  considered.  Since  the  kernels  of  the  integral 
equations  are  the  periodic  Green's  functions  and  their  normal  derivatives  on  the  surface. 
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rapid  evaluation  of  the  slowly  convergent  series  associated  with  these  functions  is  observed 
to  be  critical  for  the  feasibility  of  the  method.  The  study  has  shown  that  the  brightness 
temperature  of  the  Stokes  parameter  U  can  be  significant  in  passive  remote  sensing.  Values 
as  high  as  50  iiC  are  observed  for  certain  configurations. 

To  demonstrate  the  use  of  polarimetry  in  passive  remote  sensing  of  azimuthally 
asymmetric  features  on  a  terrain  surface,  an  experiment  was  designed  and  implemented.  A 
triangular  corrugation  pattern  was  made  on  the  sandy  soil  surface.  Polarimetric  brightness 
temperatures  are  measured  with  horizontal,  vertical,  and  45°  polarization  orientations  for 
various  observation  angles.  From  the  measured  temperatures,  absolute  values  as  high  as 
30-40  K  of  the  third  Stokes  brightness  temperatures  are  observed.  A  theoretical  analysis 
of  the  data  indicates  that  the  high  values  of  U  are  caused  by  the  azimuthal  asymmetry 
on  the  remotely  sensed  soil  surface.  It  is  also  observed  from  the  experiment  that  the 
brightness  temperatures  for  all  three  Stokes  parameters  vary  as  the  observation  direction 
varies  from  being  parallel  to  the  surface  row  structure  to  being  perpendicular  to  the  row 
structure.  The  significant  implication  of  this  experiment  is  that  the  surface  asymmetry 
can  be  detected  with  a  measurement  ofU  at  a  single  azimuthal  angle. 
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Folarimetric  Passive  Remote  Sensing  of  Periodic  Surfaces 

M.  E.  Vcytogin,  H.  A.  Yueh,  B.  T.  Shin,  and  J.  A.  Kong 

Dcputmcnt  of  Ekctrieal  Eagmcering  and  Compaict  Sdeace 
and  Rcteatch  Laboiatoty  Eleetionici 
ManachiuetU  Imtitnte  of  Techaologr 
Camlsidge,  MA  02139,  OSA 

Abstract-  The  concept  of  polarimetry  in  actire  remote  tensing  u  extended  to  passive 
remote  sensing.  The  potential  use  of  the  third  and  fourth  Stokes  parameters  V  and 
V,  which  play  an  importaat  role  in  polarimetric  active  remote  sensing,  is  demonstrated 
for  passive  remote  sensing.  It  it  shown  that,  by  the  use  of  the  reciprocity  principle,  the 
polarimetric  paimeters  of  passive  remote  sensing  can  be  obtained  through  the  solution 
of  the  assodated  direct  scattering  problem.  These  ideas  are  ^>plied  to  study  polarimetric 
passive  remote  sensing  of  periodic  surfaces.  The  solution  of  the  direct  scattering  problem 
is  obtained  by  an  integral  equation  formulation  which  involves  evaluation  of  periodic 
Green’s  functions  and  normal  derivative  of  those  on  the  surface.  Rapid  evaluation  of  the 
slowly  convergent  series  associated  srith  these  functions  is  observed  to  be  critical  for  the 
feasibility  of  the  method.  New  formulas,  whidi  are  rapidly  convergent,  are  derived  for  the 
calculation  of  these  series.  The  study  has  shown  that  the  brightness  temperature  of  the 
Stokes  parameter  V  can  be  significant  in  passive  remote  sensing.  Values  as  high  as  60  K 
are  observed  for  certain  configurations. 

J.  INTAOOUCnON 

Although  the  use  of  polarimetry  is  well  understood  in  active  remote  sensing  [1-4], 
little  attention  has  been  given  to  it  in  passive  remote  sensing.  In  active  remote 
sensing,  additional  polarimetric  information  is  obtained  from  the  third  and  fourth 
Stokes  parameters  U  and  V.  They  are  defined  as  2Re(£vJ?^}  and  2Im{£vS^) 
respectively,  where  the  brackets  denote  ensemble  avera^ng.  Therefore,  U  and 
V  are  the  measure  of  correlation  between  the  horizontal  and  vertical  components 
{Elf  and  £«)  of  the  electric  field  and  thereby  carry  further  information.  Tseng  [5] 
extended  the  concept  of  polarimetry  to  passive  remote  sensing  and  considered 
scattering  by  discrete  scatters  and  rough  sur&ces  when  the  thermal  emission  of 
third  and  fourth  Stokes  parameters  are  nonzero.  Passive  polarimetric  measure¬ 
ments  have  been  extensively  used  in  ra<fio  astronomy  to  study  planets  and  stars 
(6).  In  this  paper,  we  discuss  polarimetric  passive  remote  sensing  of  periodic  rough 
surfaces.  In  particular,  we  calculate  the  brightness  temperatures  of  U  and  V  as 
can  be  obtained  by  pasrive  measurements. 

The  reciprocity  principle  is  used  to  relate  active  remote  sensing  parameters  to 
polarimetric  pasrive  remote  sensing.  The  emisrivity  of  the  surface  for  a  particular 
polarization  can  be  obtained  as  one  minus  the  reflectivity  |7,8|.  In  order  to  be  able 
to  get  the  full  p<darimetric  information,  including  the  corresponding  brightness 
temperatures  of  V  and  V,  we  need  the  emisrivities  of  four  difierent  polarizations: 
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Classification  of  terrain  cover  using  polarimetric  radar  is  an  area  of  considerable  cur¬ 
rent  interest  and  research.  A  number  of  methods  have  been  developed  to  classify  ground 
terrain  types  from  fully  polarimetric  synthetic  aperture  radar  (SAR)  images,  and  these 
techniques  are  often  grouped  into  supervised  and  unsupervised  approaches.  Supervised 
methods,  including  both  conventional  Maximum  Likelihood  (ML)  and  more  recent  Multi¬ 
layer  Perceptron  classifiers,  have  yielded  higher  accuracy  than  unsupervised  techniques, 
but  suffer  from  the  need  for  human  interaction  to  predetermine  classes  and  training  re¬ 
gions.  In  contrast,  unsupervised  methods  determine  clsisses  automatically,  but  generally 
show  limited  ability  to  accurately  divide  terrain  into  natural  classes.  In  this  paper,  a 
new  terrain  classification  technique  is  introduced,  utilizing  unsupervised  neural  networks 
to  provide  automatic  classification,  but  employing  an  iterative  algorithm  which  overcomes 
the  poor  accuracy  of  other  unsupervised  techniques. 

Several  types  of  unsupervised  neural  networks  are  first  applied  to  the  classification  of 
SAR  images,  and  the  results  are  compared  with  those  of  more  conventional  unsupervised 
methods.  Neural  Network  approaches  include  Adaptive  Resonance  theory  (ART),  Learning 
Vector  Quantization  (LVQ),  and  Kohonen’s  self-organizing  feature  map.  Conventional 
classifiers  utilized  are  the  migrating  means  clustering  algorithm  and  the  K-means  clustering 
method.  With  both  neural  network  and  convention^  classifiers,  preprocessing  is  performed 
to  reduce  speckle  noise  and  to  stabilize  the  training  process.  Results  show  that  LVQ  is 
the  best  of  the  neural  network  techniques,  and  that  this  method  outperforms  all  of  the 
conventional  unsupervised  classifiers.  The  accuracy  of  even  the  LVQ  technique,  however, 
is  seen  to  remain  below  that  of  supervised  methods. 

To  overcome  this  poor  accuracy,  an  iterative  algorithm  is  proposed  in  which  the  SAR 
image  is  reclassified  using  a  Maximum  Likelihood  (ML)  classifier.  Training  of  the  ML 
classifier  is  performed  using  a  training  data  set  first  classified  by  the  above  unsupervised 
method,  thus,  requiring  no  human  intervention,  and  preserving  the  unsupervised  nature  of 
the  overall  classification  scheme.  The  process  is  then  repeated  iteratively,  training  a  second 
ML  classifier  using  data  classified  by  the  first.  It  is  shown  that  as  applied  to  polarimetric 
SAR  images  of  San  Francisco  and  Beaufort  sea  ice,  acquired  by  the  Jet  Propulsion  Labo¬ 
ratory,  this  algorithm  converges  rapidly,  and  significantly  improves  classification  accuracy. 
Performance  after  convergence  is  seen  to  be  comparable  to  that  obtained  with  a  supervised 
ML  classifier,  while  maintaining  the  advantages  of  an  unsupervised  technique. 
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For  polarimetric  remote  sensing,  geophysical  media  are  modeled  as  layers  containing 
randoinly  embedded  scatterers.  In  a  medium  such  as  vegetation  canopy,  the  scatterers 
can  have  various  shapes,  sizes,  and  permittivities  that  are  significantly  distinct  from  the 
background  medium.  The  model  presented  in  this  paper  will  consider  each  type  of  the 
scatterers  as  a  species  which  can  takes  on  a  shape,  size,  and  complex  permittivity  different 
from  other  species.  The  effective  permittivity  of  the  random  medium  is  derived  under  the 
strong  permittivity  fluctuation  theory  and  the  polarimetric  scattering  coefficients  are  cal¬ 
culated  for  the  layer  configuration  with  the  distorted  Born  approximation  in  the  analytical 
wave  theory  which  preserves  the  phase  information. 

The  multiple  species  in  the  random  medium  are  considered  as  randomly  oriented  ellip¬ 
soids  and  described  by  multiple  three-dimensional  ellipsoidal  correlation  functions.  The 
variances  and  correlation  lengths  of  the  correlation  functions  characterize  the  fluctuation 
strengths  and  the  physical  geometry  of  the  scatterers,  respectively.  The  strong  permittiv¬ 
ity  fluctuation  theory  is  extended  to  account  for  the  multiple  species  of  ellipsoidal  shape. 
In  the  random  medium,  a  coincidence  of  an  observation  point  with  a  source  point  gives 
rise  to  the  singularity  of  the  dyadic  Green’s  function  which  is  properly  taken  into  account 
with  exclusion  volumes  of  the  scatterers. 

Polarimetric  scattering  properties  of  a  remotely  sensed  medium  are  depicted  with  a 
covariance  matrix  whose  elements  are  polarimetric  scattering  coefficients.  The  medium 
has  a  layer  configuration:  the  top  layer  such  sis  air  is  considered  as  homogeneous,  the 
middle  layer  such  as  a  vegetation  canopy  is  random,  and  the  underlying  layer  such  as  soil 
is  a  homogeneous  half  space.  More  random  medium  layers  can  also  be  introduced  in  the 
configuration  to  account  for  weather  effect  such  as  fog  cover.  The  distorted  Born  approx¬ 
imation  is  then  used  with  the  effective  permittivity  to  obtain  the  polarimetric  covariance 
matrix.  The  result  for  the  cross-polarized  return  <7/,„  is  non-zero  even  in  the  first  order 
approximation.  Due  to  the  non-spherical  shape  and  the  random  orientation  of  the  scatter¬ 
ers,  the  correlation  coefficient  between  the  HH  and  VV  returns  has  a  magnitude  different 
from  unity  and  a  small  phase.  The  scattering  coefficients  are  also  used  to  calculate  the 
Mueller  matrix  for  synthesis  of  polarization  signatures.  The  copolarized  signature  of  the 
random  medium  has  a  rather  straight  distortion  track  and  a  recognizable  pedestal. 
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To  demonstrate  the  use  of  polarimetry  in  passive  remote  sensing  of  azimuthal  asymmet¬ 
ric  features  on  a  terrain  surface,  an  experiment  is  designed,  implemented,  and  presented 
in  this  paper. 

The  experiment  setup  is  illustrated  in  the  figure  below.  A  triangular  corrugation  pattern 
with  height  h  =  2.5  cm  and  period  p  =  5  cm  is  made  on  the  sandy  soil  surface  covering 
23  periods  by  length  t  =  160  cm  and  thickness  d  —  12.7  cm.  A  radiometer  of  15°  beam 
width  operating  at  10  GHz  is  used  in  the  experiment.  The  radiometer  is  mounted  on 
a  tripod  at  an  elevation  of  1.8  m  height  and  directed  toward  the  soil  surface  along  the 
look  direction  determined  by  azimuthal  angle  <f>  and  polar  angle  9.  Polarimetric  brightness 
temperatures  T^/j,  and  Tgp  are  measured  respectively  with  horizontal,  vertical,  and 

45°  polarization  orientations  for  9  =  20°, 30°  and  ^  from  0°  to  90°.  From  the  measured 
temperatures,  the  third  Stokes  brightness  temperatures  Ug  are  obtained. 

Absolute  values  as  high  as  30—40  K  of  the  third  Stokes  brightness  temperatures  are 
observed.  A  theoretical  aneilysis  of  the  data  indicates  that  the  high  values  of  Ug  are 
caused  by  the  azimuthal  asymmetry  on  the  remotely  sensed  soil  surface.  It  is  also  observed 
from  the  experiment  that  Tgf^  decreases,  Tg^  increases,  Tgp  decreases  to  a  minimum  at 

4>  =  45°  and  then  increases  as  <f>  takes  on  the  increased  values.  For  U g,  the  trend  is  similar 
to  that  of  Tgp.  These  general  trends  are  supported  by  our  theoretical  predictions  of  the 
polarimetric  brightness  temperatures.  The  significant  observation  from  this  experiment 
is  that  the  surface  asymmetry  can  be  detected  with  a  measurement  of  Ug  at  a  single 
azimuthal  angle. 
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Analytical  Solution  of  the  Vector  Radiative  Transfer  Equation 
with  Rough  Surface  Boundary  Condition 

H.  C.  Han,  J.  A.  Kong,  and  S.  V.  Nghiem 
Department  of  Electrical  Engineering  and  Computer  Science 
and  Research  Laboratory  of  Electronics 
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T.  Le  Toan 
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Toulouse  Cedez,  France 

Among  the  various  theoreticcJ  models  applied  to  the  study  of  electromagnetic  wave 
scatterings  from  geophysical  terrain,  such  as  snow,  ice,  and  vegetation  canopy,  the  radiative 
transfer  theory  has  drawn  intensive  attention  in  the  microwave  remote  sensing  society 
during  the  past  years.  In  most  of  the  scattering  models,  the  volume  scattering  and  the 
surface  scattering  effects  have  been  investigated  separately.  Recently,  there  has  been  a 
growing  interest  in  the  construction  of  composite  models  which  can  take  into  account  both 
types  of  scattering. 

In  this  research,  we  present  the  first  order  iterative  solution  to  the  vector  radiative 
transfer  equations  for  a  two-layer  medium  with  a  diffuse  top  boundary  and  an  irregular 
bottom  boundary  of  Gaussian  roughness.  The  Kirchhoff  approximation  and  the  geomet¬ 
rical  optics  approach  with  shadowing  correction  are  used  in  formulating  the  boundary 
conditions.  To  demonstrate  the  utilities  of  this  theory,  we  apply  our  formula  with  a  phase 
matrix  for  randomly  distributed  and  oriented  spheroidal  discrete  scatterers  to  calculate 
the  backscattering  coefficients  from  soybean  field  at  different  growing  stages  and  compare 
the  results  with  experimental  measurements.  Good  agreement  has  been  achieved  for  both 
the  co-polarized  and  the  cross-polarized  data.  It  is  observed  that  the  presence  of  the  rough 
surface  can  significantly  enhance  the  backscattering  at  small  incident  angles  and  the  levels 
of  the  cross-polarized  return.  The  polarization  signatures  calculated  based  on  the  Mueller 
matrix  shows  a  straight  distortion  track  and  an  observable  pedestal.  Numerical  comparison 
to  the  backscattering  coefficients  calculated  by  using  planar  bottom  boundary  conditions 
with  or  without  the  incoherent  addition  of  the  rough  surface  effects  are  also  made  in  this 
work. 


EXTERNAL  CALIBRATION  OF  POLARIMETRIC  RADARS  USING 
POINT  AND  DISTRIBUTED  TARGETS 

S.  H.  Yueh,  J.  A.  Kong,  and  R.  T.  Shin 

Department  of  Electrical  En^neering  and  Computer  Science 
and  Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

ABSTRACT.  Polarimetric  calibration  algorithms  using  combinations  of  point  targets 
and  reciprocal  distributed  targets  are  developed.  From  the  reciprocity  relations  of  dis¬ 
tributed  targets,  an  equivalent  point  target  response  is  derived.  Then  the  problem  of  po¬ 
larimetric  edibration  using  two  point  targets  and  one  distributed  target  reduces  to  that 
using  three  point  targets,  which  has  been  solved  before.  For  cidibration  using  one  point 
target  and  one  reciprocal  distributed  target,  two  cases  are  analyzed  with  the  point  tar¬ 
get  being  a  trihedral  reflector  or  a  polarimetric  active  radar  calibrator  (PARC).  For  both 
cases,  the  general  solutions  of  the  system  distortion  matrices  are  written  as  a  product  of 
a  particular  solution  and  a  matrix  with  one  free  parameter.  For  the  trihedral-reflector 
case,  this  free  parameter  is  determined  by  assuming  azimuthal  symmetry  for  the  dis¬ 
tributed  target.  For  the  PARC  case,  knowledge  of  one  ratio  of  two  covariance  matrix 
elements  of  the  distributed  target  is  required  to  solve  for  the  free  parameter.  Numerical 
results  are  simulated  to  demonstrate  the  usefulness  of  the  algorithms  developed  in  this 
paper. 

I.  Introduction 

In  this  paper  the  measured  scattering  matrix  corresponding  to  a  target  with  scat¬ 
tering  matrix  S  is  assumed  to  be  [1] 

X  =  cRST  (1) 

The  matrices  T  and  R  account  for  the  channel  imbalance  and  cross-polarization  coupling 
of  the  transmitting  and  receiving  ports,  respectively.  Propagation  delay  and  loss  between 
the  radar  and  the  target  is  accounted  for  by  parameter  c.  The  objective  of  polarimetric 
calibration  is  to  solve  for  the  R  and  T  matrices  using  responses  from  targets  with  known 
polarization  scattering  parameters. 

For  general  polarimetric  systems,  six  parameters  need  to  be  estimated  for  relative 
polarimetric  calibration;  three  for  each  of  the  transmitting  and  receiving  ports  including 
one  for  channel  imbalance  and  two  for  cross-talk.  To  calibrate  a  general  polarimetric 
system,  algorithms  using  three  point  targets  (in-scene  reflectors)  have  been  developed  in 
[1-5]  and  a  calibration  ^gorithm  using  one  distributed  target  and  two  point  targets  was 
presented  in  [6,7). 

For  polarimetric  radars  with  a  specific  form,  the  number  of  calibration  targets  can 
be  reduced.  For  example,  a  commonly  adopted  model  is  the  reciprocal  radar  {R  being 
the  transpose  of  T).  In  this  case,  a  esdibration  technique  using  natural  distributed  tar¬ 
gets  with  azimuthal  symmetry  and  trihedral  corner  reflectors  was  developed  18]. 

In  this  paper,  we  will  investigate  algorithms  for  calibrating  a  gener^  polarimetric 
system  using  combinations  of  point  and  distributed  targets. 

II.  Equivalent  Point  Target  Response  Using  Reciprocity 

In  this  section,  we  derive  an  equivalent  point  target  from  the  response  of  reciprocal 
distributed  targets,  utilizing  the  model  of  polarimetric  radars  defined  by  Equation  (1). 
Let  us  define  a  mapping  between  the  scattering  matrix  X  and  a  column  vector  X  as 
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where  the  superscript  t  signifies  transpose.  Likewise,  we  can  map  the  matrix  S  into  the 
column  vector  S.  Thereafter,  Equation  (1)  can  be  converted  into  a  matrix  equation 


X  =  AS  (3) 

where  A  is  the  direct  product  of  R  and  7^ 

A  =  R®T*  (4) 

For  the  case  that  5  is  the  response  from  a  distributed  target,  it  is  more  appropriate 
to  describe  the  distributed  target  by  its  covariance  matrix,  Cj.  From  Equation  (3),  the 
measured  covariance  matrices  Cx  of  Af  is  related  to  Cg  by 

Cx  =  AC,A+  (5) 

where  *+’  signifies  conjugate  transpose. 

By  assuming  the  distributed  target  is  redprocal  fSjj  =  S21)  and  its  covariance 
matrix  Cj  has  rank  3,  Klein  and  Freeman  (6,7)  derivea  the  following  equation 

G,  =  cU+"'c2,  (6) 

where  Qx  and  Qj  are,  respectively,  the  eigenvectors  of  Cx  and  Cj  with  zero  eigenvalue. 

In  addition,  S12  =  S21  means  that  columns  2  and  3  of  Cg  are  equal.  Hence, 

Q,  =  (0  -1  1  0]*  (7) 

Note  that  c'  is  an  appropriate  scaling  constant,  taking  into  account  the  arbitrary  abso¬ 
lute  magnitudes  of  eigenvectors. 

Equation  (6)  is  an  important  observation  made  by  Klein  and  Freeman  [6,7)  using 
the  reciprocity  relations  satisfied  by  the  distributed  targets.  Subsequently,  tney  carried 
out  the  matrix  inversion  and  multiplication  in  Equation  (6)  and  obtained  three  nonlinear 
equations  for  matrix  elements  of  R  and  T  [Ref.  6,  Eqs.  (19)-(21)).  Further  supplement¬ 
ing  these  three  equations  with  the  response  from  two  point  targets,  they  solved  for  the 
normalized  matrix  elements  of  R  and  T  explicitly. 

A.  Equivalent  Point  Target  Response 

In  this  paper,  we  recast  Equation  (6)  into  another  form  to  allow  for  a  general  appli¬ 
cation.  After  realizing  that  the  mapping  from  Equation  (1)  to  Equation  (3)  is  one  to  one 
in  nature,  we  can  convert  Equation  (6)  into  the  following  form 

(?+-*  =  jiRQt~^T  (8) 

where  Qg  and  Qx  are  2  by  2  matrices  corresponding  to  vectors  Qg  and  Gxi  respectively, 
through  the  mapping  defined  by  Equation  (2).  **’  denotes  complex  conjugate.  Compar¬ 
ing  Equation  (8)  with  Equation  (1),  we  can  identify  that  Q*  is  the  response  from  a 
point  target  with  scattering  matrix 

B.  Reciprocation  of  Polarimetric  Radar 

As  we  know  that  a  polarimetric  radar  is  in  general  not  reciprocal  {R  ^  T*),  there¬ 
fore  the  measured  scattering  matrix  X  is  not  symmetric  when  S  is  symmetric.  As  to  be 
shown  in  the  following,  Equation  (8)  can  be  utilized  to  symmetrize  the  measured  scatter¬ 
ing  matrix. 
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From  Equation  (8),  it  can  be  shown  that 

T  =  R*[QtQ\y'  (9) 

For  convenience,  all  the  proportional  constants  like  e'  will  not  be  written  out  explicitly 
from  now  on.  Then  substituting  Equations  (9)  into  (1),  it  is  easy  to  show  that 

Y^=XQtQ\  =  RSR*  (10) 

Note  that  if  S  is  symmetric,  so  is  Y.  After  comparing  Equations  (10)  with  (1),  we  can 
readily  interp’-^t  Y  as  the  response  observed  by  a  reciprocal  polarimetric  radar  with  T 
being  the  transpose  of  R. 

If  X  is  the  backscatter  from  a  reciprocal  distributed  target,  then  we  can  derive  from 
C'x  a  reciprocated  covariance  matrix  Cy  of  Y  which  is  given  by 


II 

(11) 

Ay  —  /  ® 

(12) 

The  results  of  this  section  show  that  the  target  reciprocity  can  be  utilized  to  derive 
an  equivalent  target  response  and  makes  a  polarimetric  radar  reciprocal. 

III.  Calibration  Using  Two  Point  Targets  and  One  Reciprocal  Distributed 
Target 

This  section  discusses  how  to  use  two  general  point  targets  and  one  reciprocal  dis¬ 
tributed  target  for  polarimetric  calibration.  Note  that  a  general  class  of  problems  where 
these  two  point  targets  are  two  reciprocal  reflectors  have  been  solved  explicitly  by  Klein 
and  Freeman  and  verified  with  the  Jet  Propulsion  Laboratory  polarimetric  sensors  [7]. 
Here,  we  will  show  how  the  equivalent  target  can  be  used  to  solve  the  problem. 

As  shown  in  the  previous  section,  we  obtain  an  equivalent  target  from  the  reci¬ 
procity.  Hence,  by  combining  this  equivalent  target  with  another  two  point  targets,  the 
problem  of  calibration  using  two  point  targets  and  one  distributed  target  is  converted  to 
that  using  three  point  targets,  which  has  uready  been  solved  [3]. 

IV.  One  Trihedral  Reflector  and  One  Reciprocal  Distributed  Target 

This  section  discusses  how  to  calibrate  a  polarimetric  radar  using  the  responses 
from  a  trihedral  reflector  and  one  distributed  target.  Note  that  van  Zyl  [8]  has  already 
provided  a  scheme  for  calibrating  a  reciprocal  polarimetric  radar  (T  =  R^)  using  one 
trihedral  reflector  and  one  distributed  target  with  azimuthal  syiiunetry.  Also  in  the  pre¬ 
vious  section  we  have  shown  that  a  general  polarimetric  radar  {R  ^  T*)  can  be  made 
symmetric  by  using  the  reciprocity,  indicating  that  van  Zyl’t  approach  can  be  directly 
applied.  However,  we  will  take  a  different  approach  in  solving  the  problem  so  that  par¬ 
tial  calibration  can  be  carried  out  without  assuming  the  azimuthal  symmetry  for  dis¬ 
tributed  targets. 

Suppose  that  the  measured  response  from  a  trihedral  is  given  by 

Xc  =  RT  (13) 

Hence, 

Yc  =  XcQtQ[  =  RR*  (14) 
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and  T  can  be  expressed  as 

T==R-^Xc  (15) 

The  general  solution  of  R  can  be  written  as  the  product  of  a  particular  solution  Rj,  and 
a  matrix  M 

=  (16) 

where  the  particular  solution  Rp  has  the  following  form: 


Rpll 

0 

■^21 

Hp22 

cos  $ 

sin^ 

—  sin  $ 

cos^ 

Note  that  6  is  allowed  to  be  complex.  After  substituting  Equation  (16J  into  Equation 
(14),  we  can  solve  for  the  following  quantities  from  the  resulting  equation 


2p21  _  l^cl2  _  ^021 

2pll  ^cll  yell 


lYe22  _  /■^2iy 
i'cll  ^RpllJ 


(19a) 

(195) 


Then  from  Equation  (15),  we  obtain  the  particular  solution  of  Tp  (=  R^^Xc)  and 
the  general  solution  of  T  (=  MTp).  Hence,  the  general  forms  of  R  and  T  have  been 
deduced  and  only  one  parameter  a  {=  tan  0)  remains  to  be  determined.  In  the  follow¬ 
ing,  we  assume  the  azimuthal  symmetry  for  the  distributed  target  to  determine  a.  At- 
imuthal  symmetry  of  the  distributed  target  implies  that 

(^$12  =  C’,13  =  C7,24  =  ^f34  =  0  (20) 

By  substituting  the  general  form  of  R  and  T  into  Equation  (1),  it  can  be  shown  that 

Xp  =  R^^XTf^  =  (21) 

Note  that  =  A/‘.  Therefore,  matrix  Xp  is  symmetric. 

Representing  Xp  by  a  vector  Xp  according  to  Equation  (2),  we  can  define  the  co- 
variance  matrix  Cp  =  E[XpX^]  which  is  related  to  (7*  and  C7»  by 


where 


Multiplying  both  sides  of  Equation  (21)  by  Af,  then  calculating  the  ensemble  av¬ 
erage  of  the  products  of  matrix  dements,  and  using  the  fact  that  Cp\2  =  Cpi3  and 
Cp34  =  Cp24,  we  obtain  the  following  equations  for  a 

0‘(Cp32  +  Ct22  -  Cp22  -  Ct32)  +  ®*(Cj,14  ~  ^tll  +  ^p44  ~  C'»4l)  =  “(^^>12  +  Cp42)  (24) 


Cp  —  ApCxAp 

(22a) 

=  AynC$A^ 

(226) 

Ap  =  R^^  9  Tp-** 

(23a) 

Am  =  9  M* 

(236) 
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“(C'psa  +  <^523  -  C'p23  -  <^j33)  +  “*(— C'pll  +  C'*14  “  Cp41  +  ^'»44)  =  ~(^pl3  +  ^'j>43)  (25) 
Note  that  a  can  be  obtained  by  utilizing  the  solution  of  the  following  equation 

aa  +  fco*  =  c  (26) 

whicli  has  the  following  solution 

“  ~  -  It.?*  provided  |c|  #  |i.|.  (27) 

The  criterion  |a|  ^  |6|  can  be  satisfied  for  a  large  class  of  distributed  targets  as  observed 
from  Equations  (24)  and  (25). 

If  we  start  with  a  polarimetric  radar  having  good  channel  isolation  (small  a),  we 
can  use  the  following  iterative  scheme  to  determine  the  value  of  a: 

Step  1.  Initialize  a  as  zero. 

Step  i.  Evaluate  Am  from  Equation  (23b)  with  the  current  value  of  a  and  calculate 
C's  =  A^CpA^^  from  Equation  (22b)  as  the  current  estimate  of  Cg. 

Step  3.  Solve  a  from  either  Equations  (24)  or  (25)  by  using  the  current  value  of  Cg. 
Step  4-  Repeat  Steps  2  and  3  until  ^i24'  and  C'3^  are  smaller  than  a  se¬ 

lected  threshold,  which  in  our  case  is  set  to  be  10~®  times 
Once  a  is  solved,  the  solution  for  R  and  T  is  complete  except  for  their  absolute 
magnitudes. 


V.  One  Polarimetric  Active  Radar  Calibrator  and  One  Reciprocal  Distributed 
Target 

This  section  discusses  how  to  use  a  PARC  and  a  distributed  target  for  calibration. 
Suppose  that  there  is  a  PARC  with  the  following  scattering  matrix  [^, 


5c  = 


0 

0 


1 

0 


(28) 


Subsequently,  substituting  Equation  (28)  into  (10)  results  in  a  equation  for  the  trans¬ 
formed  scattering  matrix  Ye-  Then  it  is  straightforward  to  show  that  the  general  solution 
of  R  has  the  following  form 

(29) 


R  =  R^M-^ 


where  the  particular  solution  is 

Rfi  =  Ri\ 

and 


1  ^cn/5cl2 
Ic22/^el2  1 


Mu  = 


1  0  ] 

0  l/0\‘ 


Here  ^  is  a  free  parameter.  In  addition,  the  particular  solution  of  T  is  obtained  from 
Equation  (8)  and  is  given  by 

Tp  =  Q^gRp^Qz^"^ 

It  can  be  shown  that  the  general  solution  of  T  can  be  written  as 

T  =  ^M-% 


(30) 

(31) 

i 

(32) 

(33) 
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Substituting  Equations  (29)  and  (33)  into  Equation  (1)  results  in 

=  R-^XT-^  (34) 

Let  Xfi  =  R~^ XT~^  and  vectorize  it  as  Xfi  by  Equation  (2).  We  can  relate  the  covari¬ 
ance  matrix  C^t  =  to  Cx  by 


—  -AfiCxX^  (35) 

where 

=  R-'  ®  T~^*  (36) 

Multiplying  Xfi  by  Rfi  from  the  left-hand  side  of  (34)  and  carrying  out  the  ensemble  av¬ 
erage  of  products  of  matrix  elements,  we  can  obtain  the  following  equations 


and 


/ ^«ii 

<^«12  _  ^«13  ^<i23 

<^j22  <^fin  C'*23  C’/il3 

^  (c^c^y 

\  <^*24  <^/i22  )  \  C1,34  ^m23  / 


0  = 


± 


Ctl4  ^p44 
Cs44  Cpl4 


(37a) 

(376) 

(37c) 


(38) 


Here  we  can  see  that  in  order  to  solve  the  parameter  /?,  it  is  necessary  to  know  one  of 
the  ratios  of  matrix  elements  of  C<  beforehand.  Also  notice  that  Equation  (37)  is  not 
applicable  to  distributed  targets  with  azimuthal  symmetry  which  will  result  in  the  case 
of  zero  divided  by  zero. 

VI.  Numerical  Simulation  and  Discussion 

In  this  section  numerical  simulation  of  algorithms  discussed  in  the  previous  two  sec¬ 
tions  will  be  presented.  The  polarization  transfer  matrices  of  the  polarimetric  radar  are 
assumed  to  be 


1  0.0426Z  -  169.5® 

0.0532Z113.6®  1.0638Z  -  86.3® 


(39) 


1  0.1042Z  -  77.8® 

0.0625Z30®  1.0417Z  -  57.9® 


(40) 


This  set  of  system  matrices  corresponds  to  the  distortion  matrices  of  the  Jet  Propulsion 
Laboratory  L-band  radar  estimated  by  a  three-point-target  approach  [9].  Note  that  the 
cross-talk  errors  are  smaller  than  —20  dB  and  the  channel  imbalances  are  around  0.5 
dB.  The  most  severe  errors  are  the  differences  in  phases  between  channels,  which  will 
significantly  distort  the  co- polarization  signature  [10]  of  trihedral  reflectors  and  make  it 
look  like  that  of  a  dihedral  reflector. 
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In  this  paper,  we  denote  channels  1  and  2  as  horizontal  (h)  and  vertical  (v)  polar¬ 
izations,  respectively.  Therefore,  the  parameters  of  a  covariance  matrix  Ct  are  defined 

as  _ 

^  's/^hvPhhhv  -^/^vhPhhvh  w^Phhvv 

Q  _  ^  '/^hvPfihhv  ^^hv^vhPhvvh  y/^hvfPhvvv 

^/^vhPhhvh  \/^f^v^vhPhvvh  y/^vhlfPvhvv 

.  Vyphhvv  V^Phvvv  y/^Plhvv  'T 

Here,  Pa0-rS  represents  the  correlation  coefRcient  of  ‘o^’  and  ‘7^*  polarizations.  For  dis¬ 
tributed  targets  with  azimuthal  symmetry,  =  0. 

Listed  in  column  1,  Table  1  are  the  covariance  parameters  of  the  distributed  target 
to  be  used  for  the  simulation  of  calibration  algorithms,  and  its  co-polarization  signature 
is  shown  in  Figure  1(a).  Due  to  the  channel  imbalance  and  cross-talk  for  the  radar  de¬ 
fined  by  Equations  (39)  and  (40),  the  observed  co-polarization  signature  would  be  that 
shown  in  Figure  l(bL  After  the  reciprocity  is  applied,  the  co-polarization  signature  cor¬ 
responding  to  the  reciprocated  covariance  matrix  Cy  is  illustrated  in  Figure  1(c),  which 
becomes  symmetric  with  respect  to  the  90-degree  prientation  angle. 

A.  One  Trihedral  Reflector  and  One  Reciprocal  Distributed  Target 

Considered  below  is  the  trihedral-reflector  approach  (see  Section  III).  Two  solutions 
were  obtained.  One  produces  the  correct  R  and  T  and  the  other  is  related  to  the  correct 
one  by 

=  r'=[j  ij]r.  (42) 

Using  these  two  solutions  to  calibrate  the  radar,  the  calibrated  covariance  parameters  of 
the  distributed  target  are  the  same  as  those  given  in  column  2  of  Table  1.  This  is  due 
to  the  fact  the  sign  difference  does  not  affect  the  responses  from  trihedral  reflectors  and 
azimuthally  symmetric  distributed  targets. 

Particular  solutions  for  the  case  considered  are  given  by 

p  _f  1  0  1  r  -f  1  0.0698Z -96.51'’’ 

■"Pl  “  [0.098/109.24°  1.06/  -  86.40°  J  ’  ~  [0.1038/22.21°  1.0482/ -  57.84° 

and 

p  _[  1  0  1  _  [  1  0.0698/ -96.51°' 

^P2  “  [0.098/109.24°  1.06/93.6°]’  0  1038/ -  157.79°  1.0482/122.16° 

Note  that  the  first  solution  is  very  close  to  the  correct  system  parameters,  in  particular, 
the  channel  imbalance.  If  the  system  is  calibrated  by  the  particular  solution,  the  cali¬ 
brated  distributed  target  parameters  are  given  in  columns  4  and  5  of  Table  1,  and  the 
corresponding  polarization  signatures  are  shown  in  Figures  1(d)  and  1(e),  respectively, 
which  are  not  visibly  distinguishable  from  the  correct  signature  shown  in  Figure  1(a). 

B.  One  Polarimeiric  Active  Radar  Calibrator  and  One  Reciprocal  Distributed  Target 

The  PARC  approach  outlined  in  Section  IV  is  considered  next.  The  particular  solu¬ 
tion  is  given  by 

p  _[  1  0.400/- 83.2“  1  r-f  1  0.1042/ -  77.8° 

“  [0.0532/113.6°  1.0  j’  ^“[0.0588/116.3°  0.9792/28.4° 

We  can  see  that  R^2i  *nd  match  correctly  with  i221  ^12  I*®®  (29)  ®nd  (33)). 

Calibration  by  the  particular  solution  alone  yields  the  parameters  of  which  are  shown 
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in  Table  2.  Note  that  aJl  the  correlation  coefficients  are  well  calibrated  except  the  phase 
of  Phhvv  Hence,  if  we  have  a  system  with  an  ideal  phase  balance  between  h  and  v  chan¬ 
nels,  then  all  the  correlation  coefficients  can  be  well  calibrated  by  the  particular  solution, 
and  T,,. 

Further  assuming  that  two  of  the  distributed  target  parameters  7  and  are 

known,  we  can  apply  Equation  (38)  for  the  calculation  of  and  to  obtdn  two  solutions. 
One  correctly  produces  the  system  parameters,  whereas  the  other  is  related  to  it  by 
Equation  (42).  By  calibrating  the  system  using  these  two  solutions,  the  calibrated  co- 
variance  matrix  parameters  are  given  in  column  2,  Table  2,  which  correctly  match  the 
original  parameters.  It  should  be  noted  that  if  correlation  coefficients,  Phhhvt  Pkhvh> 
Phvvv<  p^h-uy  are  not  zeros,  then  the  correlation  parsuneters  after  edibration  by  the 
wrong  solution  will  carry  an  opposite  sign. 

C.  Comparison  of  Trihedral  Reflector  and  PARC  Approaches 

The  results  of  the  above  two  sub-sections  indicate  that  when  the  particular  solution 
is  applied  for  calibration,  the  trihedral-reflector  approach  is  quite  robust  in  calibrating 
the  channel  imbalance  and  weak  at  detecting  the  cross-polarization  couplings.  In  con¬ 
trast,  the  PARC  approach  is  poor  at  removing  the  channel  imbalance  and  useful  in  esti¬ 
mating  the  cross-talk.  In  other  words,  the  trihedral-reflector  approach  is  useful  when  the 
system  has  relatively  small  cross-talk,  and  the  PARC  approach  is  useful  when  the  sys¬ 
tem  has  good  channel  balance,  if  no  further  assumption  about  the  distributed  targets  is 
made. 

We  should  also  note  that  the  alignment  of  point  target  is  less  critical  in  the 
trihedral-reflector  approach  than  the  PARC  approach,  since  the  PARCs  needs  to  be 
carefully  oriented  in  order  to  avoid  artificial  cross-talks.  However,  the  PARC  approach 
can  provide  a  larger  signal  to  noise  ratio  than  the  trihedral-reflector  approach. 
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Table  1.  Original,  observed,  and  calibrated  covariance  matrix  parameters  by  using 
the  responses  from  a  trihedral  reflector  and  azimuthal  symmetric  distributed  targets. 


Covariance 

Original 

Observed 

Calibrated 

Matrix 

and  Fully 

(distorted) 

by 

Parameters 

Calibrated 

and 

0.1132 


32 


0.6404 


0.6Z172.6* 


Table  2.  Original,  observed,  and  calibrated  covariance  matrix  parameters  by  using  the 
responses  from  a  PARC  and  azimuthal  symmetric  distributed  targets  with  known  7  and 
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Abstract-  A  multivariate  iC-distribution,  well  supported  by  experimental  data,  is  pro¬ 
posed  to  model  the  statistics  of  fully  polarimetric  radar  clutter  of  earth  terrain.  In  this 
approach,  correlated  polarizations  of  baclcscattered  radar  returns  are  characterized  by  a 
covariance  matrix  and  homogeneity  of  terrain  scatterers  is  characterized  by  a  parameter 
a.  As  compared  with  C-,  L-  and  P-band  polarimetric  SAR  image  simultaneously  mea¬ 
sured  by  Jet  Propulsion  Laboratory  (JPL)  on  Mt.  Shasta,  it  is  found  that  a  appears  to 
decrease  from  C-  to  P-band  for  both  the  forest  and  burned  areas. 

L  INTRODUCTION 

Polarimetric  terrain  backscatter  data  observed  with  satellite  and  airborne  sjm- 
thetic  aperture  radars  (SAR)  have  demonstrated  potential  applications  in  geologic 
mapping  and  terrain  cover  classification  [1-7].  In  previous  publications  on  this 
subject,  Gaussian  statistics  have  been  frequently  assumed  for  the  radar  return 
signals’to  build  the  Bayes  terrain  classifier  [1-3).  However,  abundant  experimen¬ 
tal  evidence  shows  that  terrain  radar  clutter  is  non-Gaussian,  i.e.,  non-Rayldgh 
in  amplitude  distribution  [8-10]. 

Among  many  non-Gaussian  statistics,  the  AT-disttibution  has  proven  to  be 
useful  in  characterizing  the  amplitude  distribution  of  electromagnetic  echoes  from 
various  objects  [9-15],  including  diverse  ground  surfaces  [9],  sea  surface  [10]  and 
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CLASSIFICATION  OF  EARTH  TERRAIN 
USING  POLARIMETRIC  SYNTHETIC 
APERTURE  RADAR  IMAGES 

J.  A.  Kong,  S.  H.  Yueh,  H.  H.  Litn, 

R.  T.  Shin,  and  J.  J.  van  Zyl 


6.1  Introduction 

6.2  Supervised  and  Unsupervised  Classifications 

a.  Supervised  Classification  Procedure 

b.  Unsupervised  Classification  Procedure 

c.  Classification  Error  Analysis 

d.  Sununaiy 

6.S  Maximum  Contrast  Enhancement 

a.  Polarimetric  Matched  Filter  and  Contrast  Ratio 

b.  Optimal  Polarimetric  Matched  Filter  Required  to  Obtain 
Maximum  Contrast  between  Two  Scattering  Classes 

c.  Optimal  Recdving  Polarization  State  for  a  Fixed 
l^ansmitting  Polarization 

d.  Results  and  Discussion 
Acknowledgments 
References 


0.1  Introduction 

Classification  of  earth  terrain  within  an  image  is  one  of  the  many 
iiiq>ortant  iq>plications  of  polarimetric  data.  A  systematic  classification 
procedure  will  place  the  classification  process  on  a  more  quantitative 
kvd  and  reduce  the  amount  of  photo-interpretation  necessary  [Evans, 
1986].  Single  feature  and  mult^quency  classifications  have  been  used 
ht  the  past,  but  dasdfication  can  now  be  applied  to  fully  polarimet¬ 
ric  data  which  have  become  available  due  to  recent  devdopments  in 
radar  technology.  It  has  been  shown  that  Bayes  dassification  using 
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Note  tliat  dusification  performaace  is  independent  of  the  nnrmAli^i,. 
tion  scheme  since  the  normalisation  factor  does  not  appear  in  the  above 
expression. 
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Abstract-  Strong  permittivity  fluctuation  theory  if  tifed  to  solve  the  problem  of  scatter' 
ing  from  a  medium  composed  of  completely  randomly  oriented  scatterers  under  the  low 
frequetrey  limit.  Based  on  Finkel’berg’s  approach  (2|3]i  Gaussian  statistics  is  not  assumed 
for  the  renotmalixed  scattering  sources.  The  effective  permittivity  is  obtained  under  the 
low  frequency  limit  and  the  result  is  shown  to  be  isotropic  due  to  no  preferred  direction 
;  in  the  orientation  of  the  scatterers.  Numerical  results  of  the  effective  permittivity  are 
illustrated  for  oblate  and  prolate  spheroidal  scatterers  and  compared  with  the  results  for 
spherical  scatterers.  The  results  derived  are  shown  to  be  consistent  with  the  discrete  scat- 
'  terer  theory.  The  effective  permittivity  of  random  medium  embedded  with  nonspherical 
'  scatterers  shows  a  higher  imaginary  part  than  that  of  spherical  scatterer  case  with  equal 
correlation  volume.  Under  the  distorted  Born  approximation,  the  polarimetric  covariance 
matrix  for  the  backscattered  electric  field  it  calculated  for  the  half-space  randomly  ori¬ 
ented  scatterers.  The  nonspherical  geometry  of  the  scatterers  shows  significant  effects  on 
the  cross-polarized  backscattering  returns  <r»,  and  the  correlation  coefficient  p  between 
HH  and  VV  returns.  The  polarimetric  backscattering  scattering  coefficients  can  provide 
useful  information  in  distinguishing  the  geometry  of  scatterers. 


I.  INTRODUCTION 

There  has  been  a  continually  growing  interest  in  scattering  from  random  media 
embedded  with  nonspherical  particles  because  the  particles  in  rain,  ice  crystals, 
fog,  snow,  leaves,  etc.,  are  nonspherical.  The  radiative  transfer  theory  has  been 
used  to  study  a  layer  of  randomly  positioned  and  oriented  nonspherical  particles 
overlying  a  homogeneous  dielectric  half  space  [17,15].  A  second  approach  based  on 
wave  theory  and  Foldy’s  approximation  [6] ,  which  is  valid  only  for  sparse  concen¬ 
trations  of  scatterers  [5]  has  also  been  developed.  For  nonsparse  concentrations  of 
i  scatterers,  incorporating  the  correlation  of  scatterer  positions,  the  quasicrystalline 
approximation  and  distorted  Born  approximation  have  been  used  to  calculate  the 
i  mean  field  and  the  incoherent  field  [12,14). 

I  For  a  medium  with  large  spatial  permittivity  variations,  such  as  water  mois- 
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ture  in  vegetation  and  brine  inclusions  sea  ice,  strong  fluctuation  theory  has  been 
developed  for  uniformly  aligned  spheroidal  scatterers  [11,13,16].  Because  the  sin¬ 
gularity  of  the  dyadic  Green’s  function  is  properly  taken  into  account  in  the  renor¬ 
malization  method,  the  theory  is  applicable  to  both  small  and  large  variances  of 
permittivity  function.  However,  the  orientations  of  scatterers  for  natural  terrain 
media  are  in  general  random,  and,  in  this  paper,  we  extend  the  theory  to  take 
into  account  the  effect  of  random  orientations. 

In  the  derivation  of  strong  fluctuation  theory,  a  Gaussian  random  process  is 
ususJly  assumed  for  the  renormalized  scattering  sources  [9,10,11,16].  The  Gaus¬ 
sian  statistics  assumption  is  employed  to  factor  the  higher  order  moments  of  the 
fluctuations  in  terms  of  the  products  of  second  moments,  and  a  Feynman  diagram¬ 
matic  representation  of  Dyson’s  equation  for  the  mean  field  is  derived.  Called  the 
bilocal  approximation,  the  mass  operator  in  the  form  of  an  infinite  series  is  ap¬ 
proximated  by  the  first  term.  The  effective  permittivity  is  derived  through  the 
constitutive  relation  between  the  mean  electric  field  and  the  mean  dielectric  dis¬ 
placement.  In  the  case  of  non-Gaussian  fluctuations,  Finkel’berg  [2,3,4]  introduced 
the  one-group  approximation  for  the  mass  operator.  When  the  fluctuations  are 
Gaussian,  the  one-group  approximation  reduces  to  the  Gaussian  case.  In  this 
paper,  Finkel’berg’s  approach  is  generalized  to  randomly  oriented  nonspherical 
scatterers.  The  polarizability  operator  is  also  in  the  form  of  an  infinite  series. 
The  first  term  of  the  series  is  the  same  for  either  Gaussian  or  non-Gaussian  fluc¬ 
tuations.  This  indicates  that  if  the  series  is  approximated  by  the  first  term,  the 
bilocal  approximation  for  the  effective  permittivity  gives  rise  to  the  same  result 
for  Gaussian  and  non-Gaussian  fluctuations. 

In  Section  11,  the  effective  permittivity  for  an  unbounded  medium  embedded 
with  isotropically  oriented  nonspherical  scatterers  is  derived  using  the  strong  fluc¬ 
tuation  theory.  The  effective  permittivity  tensor  is  expressed  in  terms  of  an  infinite 
series.  The  mean  field  is  shown  to  propagate  inside  a  medium  with  the  effective 
permittivity.  The  series  is  truncated  to  the  lowest-order  correction  term  under 
the  bilocal  approximation.  An  analytical  expression  for  the  effective  permittivity 
is  given  for  a  specific  correlation  function  under  the  low  frequency  approximation. 
In  Section  III,  numerical  results  for  the  effective  permittivities  for  spheroidal  scat¬ 
terers  are  compared  with  the  discrete  scatterer  theory  in  the  common  regime  of 
validity.  In  Section  IV,  under  the  distorted  Born  approximation,  the  scattering 
from  a  half-space  with  ellipsoidal  scatterers  is  solved.  The  covariance  matrix  of 
the  polarimetric  backscattering  coefficients  is  obtained.  Because  of  the  azimuthal 
symmetry  of  the  orientation,  four  out  of  nine  elements  of  the  covariance  matrix 
are  shown  to  be  zero  [1].  In  Section  V,  the  polarimetric  backscattering  coefficients 
are  illustrated  for  spheroidal  scatterers  with  different  axial  lengths. 

n.  STRONG  FLUCTUATION  THEORY 

Consider  an  unbounded  space  filled  with  randomly  distributed  and  oriented  scat¬ 
terers.  Assuming  the  medium  is  characterized  by  the  permittivity  ((r),  the  wave 
equation  for  a  time  harmonic  electric  field  1?(F)  is 
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l2  ^(’’)  I?/— \ 


y  X  V  X  £;(r)  -  fcg  ^  £(?)  =  0 
«o 


where 

(ft,  re  scatterer 

^  ^  I  «i>,  re  background  ^  ^ 

Introducing  tg  which  will  turn  out  to  be  the  lowest  order  term  of  the  effective 
permittivity  by  solving  (18),  we  can  recast  (1)  into 

V  xVx£(r)~t^'-i£(r)  =  k^^^^^^E(r)  (3) 

«0  «0 

By  introducing  the  dyadic  Green’s  function 

y  X  y  X  5j(f,r')  ^g{f,r)  =  lS{f  -  r")  (4) 

*0 

(3)  is  transformed  into  the  foUowing  integral  equation; 

E(f)  =  Eo  +  kl  I  dr'dg(f,f') .  •  ECf')  (5) 

Eq  is  the  applied  electric  field.  In  consideration  of  the  singularity  at  the  source 
point,  the  dyadic  Green’s  function  is  decomposed  into  two  parts; 

5,(r,f’')  =  FS^,(F,r')-%i«(r-f')  (6) 

*0 

where  PS  denotes  the  principal  value  and  5 (F)  takes  into  account  the  shape  of 
exclusion  volume  at  r  (16]. 

By  substituting  (6)  into  (5),  the  following  equation  is  obtained; 

f  (r)  =  Eo  +  kl  I  df'PSdg(f,r') .  f(r')  •  F(f‘)  (7) 

where  F(r)  is  the  external  field  and  €o((r)  ■  F{f)  corresponds  to  the  induced 
polarization  density: 

F(r)  =  [j  +  f(r)  m  (8a) 

Expand  the  solution  of  the  integral  equation  (7)  in  terms  of  infinite  series, 

T{f)  =£o  +  ''o  J •^tPS'Gg(f,fi)  1(71)  ■’Eoifi) 

+  ^0  / J  •tr2PS 5y(F,ri)  •  f(ri)  ■  PS  ^g{ri,72)  •  f(Fj)  •  So(r2) 
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+  fc*  J  dri  jdr2  J  drsPS  5^,(r,ri)  •  |(ri)  •  PS  Gg{ri,r2)  ■  ((^2) 


■PSGg(r2,r3).((f3)Eo(r3) 


+  ■■■  (9) 

The  above  equation  can  be  represented  in  an  operator  form 

F(r)  =  £0+  PSWg(r,r3)-((r3)Eo(n) 

+  PS  Gg(r,fi)  .  ((rj)  •  PS  Gg(ri,r2)  ■  ((fj) 

+  PS  Gg(r,ri)  ■  f(ri)  •  PS  •  T(»"2)  ’  Gg(r2,r3)  ■  ^(fs)  •  Eoirs) 

+  •••  (10) 

Taking  the  ensemble  average  of  (10),  we  obtain 

(F)  =  (7  +  (PSS,  •!)  +  (PSWg-t  PSUg  ■!) 

+  (PSWff-hPS^g  (  PS^g  i)+  -)  -Eo  (11) 

Multiplying  (10)  by  ((7)  and  taking  the  ensemble  average, 

{t-T)=({t)  +  {i-PSdg-()  +  {tPSdg.j.PSdg  t) 

+  (f  •  J>5 -f  ■  P5 S',  •?)  +  •••)  -  Eo  (12) 

Through  cancelling  ISq  from  (12)  by  using  (11),  the  following  equation  is  obtained: 

«o(!-y)-*o!«//(7)  (13) 

where 


=  (^>  +  ((^  •  -  (^>  •  ■  ^)) 

+  ((!  •  P5S,  •  pss,  •!)  -  (!  •  psdg  ■ !)  •  {psGg  ■ !) 

-  (f)  -(PSUg-f-  PSUg  ■  !)  +  (f)  •  (PSWg  ■()-(PSWg-  ()) 

+  -••  (14) 

Note  that  (13)  relates  the  mean  induced  polarisation  density  to  the  mean  exter¬ 
nal  field.  Making  use  of  (8),  we  can  relate  the  mean  induced  polarization  density 
and  the  mean  external  field  to  the  mean  dielectric  displacement  (S)  and  the  mean 
electric  field  ("E)  .  Thus  the  following  constitutive  relation  can  be  obtained: 

<2?(i')>=f,/,(5(f))  (15) 

?€//  =  «,7  -t-  CO  (?  -  ft//  (f )  )  (t// 


where 


(16) 
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is  the  effective  permittivity  tensor  of  the  random  medium.  Note  that  in  the 
derivation  of  (15)  and  (16),  the  fact  that  the  D  and  E  are  independent  of  S  or 
the  shape  of  exclusion  volume  has  been  used  [18]. 

By  taking  the  ensemble  average  of  (3)  and  making  use  of  (15),  the  meem  field 
can  be  shown  to  satisfy  the  following  wave  equation: 

VxVx(E)-kl‘-^(E}  =  0  (17) 

*0 

This  indicates  that  the  mean  field  inside  the  random  medium  propagates  in  a 
homogeneous  medium  with  the  effective  permittivity  tensor  ?*//• 


Effective  Permittivity  Under  Low  Frequency  Approximation 

In  evaluating  the  effective  permittivity  (16),  eg  and  need  to  be  computed. 
Note  that  gives  the  lowest  order  correction  to  the  effective  permittivity  as 
indicated  by  (l4)  and  (16).  In  this  paper,  tg  will  be  adjusted  so  that  the  first 
term  of  {«//,  (14),  is  zero,  i.e., 

(!)=0  (18) 

Also,  the  infinite  series  (14)  will  be  approximated  by  the  lowest-order  term 

Iff  ^  {t  PS^g -t)  (19) 

The  above  approximations  correspond  to  the  bilocal  approximation  for  strong 

fluctuation  theory  (17). 

Assuming 

|(!e//)i;l  <  1  (20) 


the  effective  permittivity  of  the  medium  under  the  low  frequency  approximation 
is  given  by  [Appendix  A] 

=  =(0) 

W/ - + '0  W/  (2^) 

=  {0) 

where  is  given  in  {AG).  Explicitly, 

=  jd9d<l,P{e{r),4>{f)) 

(22) 

(^0^  J  !(?)  •  Gg{r,ff)  .  f(r')  d/  +  ((F)  .  5(r)  •  ((f)|0(r),  ^(r)^ 

where  8{f)  and  ^(r)  represent  the  corresponding  orientation  angle  of  the  scat- 
terer  at  r  and  P{6,  ^)  corresponds  to  the  probability  density  function  of  the 
orientation. 

Let 

(ff(n]  M  <>(r),^(r)}=r^m>:M  (’'.5^)  (23) 

\l  Jill  / 
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and  assume  the  normalized  CMrelation  function  for  the  cross-correlation  between 
the  components  of  ^(f)  and  ((r^)  to  be 


V 


(24) 


(25) 


^  iimj;  (^»  ^) 

Using  Einstein’s  convention 

^^0  y  •  ^gir,r)  ■j{f‘)df'  +  f  (r)  -  ?(r)  •  f(f)|«(f ),  ^(r) 

=  j  ^  (r, r)  y (^>’^)  \p9  [5],  I 

5  will  be  chosen  in  such  a  way  that  jf  approaches  zero  when  the  frequency 
goes  to  zero. 

y  (26) 

In  this  manner  tg  corresponds  to  the  effective  permittivity  at  very  low  frequency. 
Assuming  the  scatterers  are  of  spheroidal  shape,  let 


^0  y -^=0,^=0  = 


and 


/»  0  0 

0  /p  0 

0  0/, 


=  ,  5p  0  0 

=  0  S.  0  I 

'0  (T  5  J 


=  —  lim 
«— *0 


/p  0  0 

I  0  Ip  0  I 

Lo  6  h\ 


Defining  the  local  coordinates  with  the  following  transformation  [Fig.  1) 
«3  =  *  sin  9  cos  ^  +  y  sin  ^  sin  ^  -t-  z  vOS  9 
*1  =  *  sin  ^  —  y  COS  (ft 

*2  =  *  cos  9  cos  ^  -f  y  cos  9  sin  ^  —  z  sin  5 
such  that  ej  is  in  the  direction  of  the  axis  of  the  scatterer.  Therefore, 

f  =  f  f,  .  f 
5  =  .  f 

where  the  effective  scattering  source  in  the  scatterer  coordinates  is 


(27) 

(28) 

(29a) 

(296) 

(29c) 

(30o) 

(306) 


<  -  <0 


<0  +  5V(«  -  tg) 
0 


e  -  c. 


*0  +  ~  *y) 

0 


0 

J-ZlL. 


<0  +  5,(«  -  €g)  J 


(31) 
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and  the  coordinate  transformation  matrix  is  given  by 

sin  4>  —  cos  4>  0 

2^  _  I  cos  0  cos  ^  cos  d  sin  ^  —  sind 

sin  d  cos  ^  sin  0  sin  ^  cosB 


z 


Figure  1.  Coordinate  system  and  scatterer  geometry. 
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(32) 


Because  the  scatterers  are  assumed  to  be  isotropically  oriented,  the  probability 
density  function  of  the  orientation  is 

=  (33) 

After  straightforward  algebraic  manipulation,  the  following  result  is  obtained: 

(|(r))  =  J 

=  [/  (I  („ + i  <,.)  +  (1  -  /)  (I  a, + i  a,)]  ?  (34) 

and 
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(i//  =  {/  [I  (Ip  +  5p)  +  ^  (I-  +  s,) 


+(i-/) 


5e6p(^P  +  ^p)+^f6z(^^  +  5x) 


where 


t«  -  <9 


(35) 

(зба) 

(збб) 
(36c) 
(36d) 

and  /  is  the  fractional  volume  of  the  scatterers.  The  subscripts  a  and  6  represent 
scatterer  and  background,  respectively. 

For  a  given  correlation  function  (24),  e.g.,  (37),  Sp  and  St  ait  related  to  tg 
by  (28).  The  equation  for  tg  is  obtained  by  letting  (34)  equal  to  zero  [see  (18)] 
for  given  Sp  and  St-  These  three  coupled  equations  constitute  the  basis  for  the 
solution  of  tg  1  ^p>  and  Sf 

In  this  paper  the  correlation  function  in  the  local  coordinates  (ei,e2,e3)  is 
chosen  to  be 


«0 

+  sp(it  ~ 

‘ff) 

—  £j 

«0 

+  ■5,(e,  — 

fb-fg 

«0 

+ 

1 

^b-tg 

«0 

+  ■^*(*6  - 

■R«=0.  ^=o(»')  =  «P 


(37) 


For  this  correlation  function,  analytical  expressions  of  Sp,  St,  Ip,  and  It  are 
given  in  Appendix  B. 


m.  NUMERICAL  RESULTS  OF  EFFECTIVE  PERMITTIVITY 


In  this  section,  we  compare  the  effective  permittivities  obtained  using  strong  fluc¬ 
tuation  theory  (21)  and  discrete  scatterer  theory  (Cl)  for  randomly  oriented 
spheroids  under  the  low  frequency  and  small  fractional  volume  limits.  The  effec¬ 
tive  permittivity  obtained  with  the  discrete  scatterer  theory  is  given  in  Appendix 
C.  The  effective  permittivities  are  plotted  as  a  function  of  frequency.  The  parame¬ 
ters  of  the  correlation  function  (37)  are  related  to  the  discrete  scatterer  parameters 
in  the  following  manner.  The  correlation  volume  is  selected  to  be  the  volume  of 
the  discrete  scatterers,  i.e.. 


4w 

T 


61*1, 


(38) 


and  the  equi-correlation  surface  is  chosen  according  to  the  shape  of  the  scatterer 


Ip 

h 


a 

c 


(39) 


Imaf»n#ry  pari  of  r, 
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Figure  2.  Imaginary  pari  of  as  a  function  of  frequency,  (a)  Ip  = 

linni  and  /.  =  0.0G25mm,  (b)  Ip  =  2mm  and  /.  =  0.25  mm, 
(c)  Ip  -  1mm  and  I.  =  1mm,  (d)  Ip  =  0.5mm  and  i;  = 
4  mm,  and  (c)  Ip  =  0.25  mm  and  =  16  mm.  The  dashed 

line  corresponds  to  discrete  scalterer  theory.  The  solid  curve 
corresponds  to  strong  fluctuation  theory. 
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Figure  2.  Continued. 


rrasuaner  <OHa)  Fraquancx  (OHa) 


Figure  3.  Real  part  of  (gff  as  a  function  of  frequency,  (a)  Ip  =  4  mm  and 
Iz  =  0.0625  mm,  (b)  Ip  =  2  mm  and  U  =  0.25  mm,  (c)  Ip  = 
1mm  and  Iz  =  1mm,  (d)  Ip  =  0.5  mm  and  Iz  =  4  mm,  and  (e) 
Ip  =  0.25  mm  and  Iz  =  16  mm.  The  dashed  line  corresponds  to 
discrete  scatterer  theory.  The  solid  curve  corresponds  to  strong 
fluctuation  theory. 
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®<r>e«>o«o«»e«o  ^0OMe«»e«eeo 

rr«q««ac7  <0H«)  rr^utne^  (OH*) 

Figure  3.  Continued. 

The  imaginary  and  real  parts  of  the  effective  permittivities  are  shown  in  Figs.  2 
and  3,  respectively.  The  permittivities  of  the  scatterers  and  background  are  as¬ 
sumed  to  be  real  so  that  the  scattering  will  be  the  dominant  effect  compared  to  the 
absorption  effect.  The  fractional  volume  of  the  scatterers  is  0.5%.  The  correlation 
volume  (38)  is  kept  a  constant  for  all  cates.  Figures  2(a)  and  2(b)  correspond  to 
disk-like  scatterers  with  correlation  lengths  {Ip  =  4.0  mm  and  =  0.0625  mm) 
and  {Ip  =  2.0  mm  and  Iz  =  0.25  mm),  respectively.  Figures  2(d)  and  2(e)  are  for 
needle-like  scatterers  with  correlation  lengths  {Ip  =  0.5  mm  and  Iz  =■  4.0  mm), 
and  {Ip  =  0.25  mm  and  Ij  =  16  mm),  respectively.  Figure  2(c)  it  for  spherical 
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scatteiers  (Ip  =  1.0  mm  and  I,  =  1.0  mm).  The  corresponding  cases  for  the 
real  parts  of  are  shown  in  Figs.  3(a)  through  3(e).  Note  that  the  real  and 
imaginary  parts  of  for  the  nonspherical  scatterers  are  larger  than  the  spher¬ 
ical  scatterer  case.  This  indicates  that  the  nonspherical  scatterer  has  a  higher 
scattering  effect  than  the  spherical  scatterer.  Also,  the  effective  permittivities 
obtained  using  strong  fluctuation  theory  and  discrete  scatterer  theory  are  in  good 
agreement  with  each  other.  The  noticeable  discrepancies  in  Fig.  2(e)  for  frequen¬ 
cies  greater  than  3  GHz  (which  msikes  kit  >  1)  is  due  to  the  fact  that  the  low 
frequency  approximation  is  used  to  derive  results  for  both  approaches. 

IV.  backscattering  from  a  half  space  random  medium  ^ 


Considering  an  incident  electric  field  impinging  on  a  half  space  filled  with  ran¬ 
domly  distributed  and  oriented  ellipsoid  scatterers  [Fig.  4]: 

®0«  =  [M-*:0si)^W  +  t»(-*:02f)-Bvi]  exp  (ifcp,  -  p  -  iktiz)  (40) 

The  unperturbed  field  inside  the  medium  is 

=  [r*(l:2.)A, (-!:,„)£*.■  +  T,(kti)vi(-kui)Eti]  exp  (ikpi  ■  p  -  ik^iz) 


n(kt)  = 


itx  +  ifci, 

r  rt  t  2<il!2  *0 

€ikt  +  €okit  ki 

The  scattered  field  in  the  far-ffeld  is  ^ven  as 


E‘=klf  <tpi  f  dziU(,i{r,ri).((fi)-F(fi) 

•/-oo  y— oo 

where  Uoi  is  the  dyadic  Green’s  function  given  by  |16j 

^oi(ry)  =  ^  Uoi(kt)h(kt)hi(kit)  +  ^yoj(kt)i>(kt)vi(k,t) 


where 


exp(-iki  rj) 

(44) 

V  2kt 

(45a) 

y  2eikt 

(456) 

kp 

(46a) 

-xkgkt  -  iikykt  +  ikj 

M**)  *  r  t 

(466) 
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*p 

(46c) 

.  f.  .  -^kxku  -  ykyku  +  ikl 

n(kiz)  = - - - 

k\kp 

(46d) 

kx  =  to  cos  sin  9, 

(47a) 

ky  =  ko  sin  <i>,  sin  9, 

(476) 

kt  =  fco  cos  9, 

(47c) 

kzi  —  ko  cos  sin  9^ 

(47d) 

kyi  =  ko  sin  sin  Pj 

(47e) 

kzi  —  ko  cos  9i 

(47/) 

Z 


Figure  4.  Configuration  for  scattering  from  half-space  with  randomly  ori¬ 
ented  spheroids. 

Applying  the  distorted  Born  approximation  to  crdculating  the  backscattering 
coefficients,  we  let 


F(f)  ~  £W{r)  (48) 

Therefore,  by  substituting  (48)  into  (43),  the  scattered  field  is  given  as 

^  =  '^’^0  «P  (-‘*1  •  ^l)  «P  •  Pi  -  »*lst*l) 

^XoiM**)  Al(*:ii)  +  ^0iv(*x)  Vi(*1j)| 

•  !(n)  •  lTHhi{-kui)E^i  +  r,ti, (-*,„)  £,.]  (49) 
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The  scattered  field  is  then  expressed  in  the  following  form: 


where  and  Ev,  are  the  horizontally  and  vertically  polarized  electric  fields 
respectively. 

Note  that  for  backscattering  from  a  reciprocal  medium  ff^^  equals  to  Th< 
polarimetric  covariance  matrix  in  the  backscattering  direction  is  denoted  as  [1] 

^hh  ^hhhv  ^hhvv 
^'hhkv  *^hvvv 

^kkvv  ‘'hvvv 

’{\fhhf)  (fhhfH,)  Uhhf:,)' 

= 

(/vv/w)  ihM  (l/.vl') 

Let  the  spectral  density  function  be 

=  J  (52 

The  backscattering  scattering  coefficients  are  given  by 
<rKh  =  I^OiTfcl^  1^4,  J’dO  ~ 

|/|{«p  +  ({»*  -  M  sin*  9  sin*(^j  -  ^)| 

+  (1  -  /)  +  {(bz  -  (bp)  sin*  9  sin*{^<  -  .^)|*|  'I'  (2kpi,9, 4>)  (53c 

/L  ,  k  .  \2]2 

/  (»p  -r  (O*  -  ($p)  si”  ^  cos(^i  ~  ^  j 


tr*,=irkS|jro,r,|*jr’'d^jf' 
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sin*  6  sin*(^i  —  sin  5  cos(^{  —  ^)  + 

*1 


and  the  cross  correlation  between  polarizations 

r2»  /•» 


^cos^l  Y 


i2k^,9,4>) 

(53c) 


<'khvv  =  ^k*  XoiTf,  Yo\T:  I  ’d4>£de  ^ 

I  /  \itp  +  {(iz  -  (zp)  sin*  tf  sin*{^j  -  ^)J 

^  (zp  +  -  izp)  cos(^i  ~  cosS^ 

+  (1  -  /)  [ifcp  +  (4*  -  (bp)  s‘n^  ^  sin*(.^,  -  ^)j 

^  (bp  +  ((bx  -  (bp)  sintf  cos(<ij  ~  4>)  +  cos  8^ 


^{2kpiA<(>) 

(54a) 

•^hhhv  ~  0 

(546) 

^Kvvv  ~  ® 

(54c) 

where 

1 

{kz  +  2kizi){kz  - 

^  {2k^i,e,<t>)  = 

r  dkz 

/ — oo 

(- 2jfcJ  i^^pitkx) 

(55) 

Under  the  low  frequency  approximation, 

i{2k^ 

(56) 

Thus 

<'hh  = 

:n7^l^o.7Xl*  / 

2*u.  V 

{|{«p[  +  g|((»*  -  (xp)l^ 

+  jR*  [(jp((><  —  (zp)\  j 

r  -  Mr  +  5R' 

[«6*p(ftx  -  (fcp)]  }  )  *9.^(0) 

jr*iki  , 

<r„  =  ^  lybiTVi* 


]^j  +  5  l({«  ~  {»/>)!*  +  [CpU**  ~  (zp)] 


”  i<«  -<*"1*} + +  gi«»x  -  Ml’ 
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+  n  ~  *«,^{0)  (57fc) 

~  2k^  '-*0l^‘'l  ^  l^-' |(6x  -  4»p)| 

+  (l-/)|«6x-ap)|"]*M(0)  (57c) 

hH.v  =  ^  JToiTfcKo'ir;  (^y  ^/||6p|'  +  -  M 

+if,.(f„-u-[(^)%(^)]' 

+  ilf„-<„p[(^)  +(^)  ]  } 

+  (1  -  /)|  I^6p|*  +  j(bp  ((iz  -  (bp)  +  jibp  ((bz  -  (bp)*  (j^)^  + 

+  ^Ifix  -  (bp\^  [(^)  +  (^)  j  I)  ♦m(O)  (57«f) 


V.  NUMERICAL  EXAMPLES  OF  POLARIMETRIC  COVARIANCE  MATRICES 

In  this  section  the  results  of  Section  IV  are  illustrated  using  the  parameters  ap¬ 
plicable  to  microwave  remote  sensing  of  vegetation  and  forests.  The  orientations 
of  nonspherical  scatterers  are  assumed  to  be  completely  random.  The  correlation 
function  (37)  is  used  to  evaluate  all  of  the  numerical  results.  Because  four  elements 
of  the  covariance  matrix  are  zero,  it  is  more  convenient  to  write  the  polarimetric 
covariance  matrix  (51)  into  the  following  form  [1] 

_  10  Pyff' 

^  =  <thh  0  '  0  (58) 

p*v/7  0  7 

where  <  and  7  correspond  to  the  ratios  Vhvl^hh  <^vvl<^hh>  respectively,  and 
p  is  the  complex  correlation  coefficient  of  of  HH  and  VV  polarization  returns. 
The  polarimetric  covariance  matrix  contains  very  important  information  in  the 
classification  of  various  terrain  elements  [7]. 

The  effects  of  nonspherical  geometry  on  polarimetric  covariance  matrices  are 
shown  in  Tables  1  and  2.  Included  are  results  for  five  types  of  scatterers  discussed 
in  Section  III.  Note  that  for  the  spherictJ  scatterer  case  is  bigger  than  the 
nonspherical  scatterer  cases.  This  it  due  to  the  fact  that  the  scattering  loss  (imag¬ 
inary  part  of  <«//)  for  the  spherical  scatterer  case  it  the  smallest  among  all  cases. 
The  wave  can  penetrate  deeper  into  the  random  medium  with  smaller  scattering 
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loss.  As  a  result,  for  a  half-space  configuration,  the  total  scattering  increases.  Also 
note  that  the  ratio  <rvvl<^hh  close  to  1  for  all  cases.  This  is  because  all 

the  scatterers  are  completely  randomly  oriented  so  that  the  number  of  scatterers 
which  contribute  to  the  HH  return  is  the  same  as  the  number  of  scatterers  for 
the  VV  return.  The  reason  why  the  ratio  is  slightly  greater  than  one  is  because 
the  transmission  coefficients  for  vertical  polarization  (42b)  and  (456)  are  slightly 
greater  than  those  of  horizontal  polarization  (42a)  and  (45a). 


Shape 

Ip  (mm) 

Ig  (mm) 

<'hh  (48) 

"vv/irtth 

P 

oblate 

4 

0.0625 

-24.6 

1.002 

0.108 

0.784 

oblate 

2 

0.25 

-24.2 

1.0006 

0.086 

0.829 

sphere 

1 

1 

-23.5 

1.0 

o 

d 

1.0 

prolate 

0.5 

4 

-25.1 

1.0001 

0.214 

0.572 

prolate 

0.25 

16 

1 

-26.0 

1.0002 

0.265 

0.470 

Table  1.  Covariance  matrix  elements  for  five  types  of  scatterers.  Fre¬ 
quency  =  1.1  GHz,^i  =  40°,  /  =  0.005,  tg  =  (30.8  -f  tl.8)eo> 
and  <(  =  CQ. 


Shape 

If  (mm) 

Iz  (mm) 

ohh  (4B) 

^vv/^Ui 

‘'hv/<’hh 

f 

oblate 

4 

0.0625 

-24.8 

1.006 

0.107 

0.787 

oblate 

2 

0.25 

-24.4 

1.002 

0.0»' 

0.830 

sphere 

1 

I 

-23.5 

1.0003 

1.0 

prolate 

0.5 

4 

-25.3 

1.001 

0.212 

0.577 

prolate 

0.25 

16 

-26.1 

1.002 

0.262 

0.477 

Table  2.  Covariance  matrix  elements  for  five  types  of  scatterers.  Fte- 
quency  =  1.1  GHz, =  40®,/  =  0.01,  c«  =  (30.8 -1- tl.8) <o, 
and  =  <0- 

The  nonspherical  geometry,  however,  has  significant  effects  on  ^hv/^hK  P- 
The  imaginary  part  of  p  is  not  shown  in  the  table,  since  they  are  very  close  to  zero. 
It  can  be  seen  that  when  the  scatterers  become  less  spherical,  ffhvl^^hh  increases 
from  zero  and  p  decreases  from  1.  When  the  eccentricity  of  the  scatterers  gets 
higher,  the  returns  for  HH  and  VV  will  be  dominated  more  and  more  by  scatterers 
with  different  orientations  (i.e.,  different  location).  As  a  result,  the  correlation 
between  HH  and  VV  decreases.  As  a  matter  of  fact,  for  completely  randomly 
oriented  dipoles  both  the  ratio  o-hv/'^hh  '^nd  p  are  equal  to  1/3.  This  is  consistent 
with  the  trend  shown  in  Table  1. 

The  results  shown  in  Tables  1  and  2  are  numerically  integrated  from  (53)  and 
(54).  Because  the  integrands  are  relatively  smooth  analytic  functions,  Simpson’s 
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rule  is  used  for  numerical  integration.  The  total  computation  time  for  one  co- 
variance  matrix  is  about  four  seconds  on  our  VAXstation  3500.  The  approximate 
formula  given  in  (57)  are  also  evaluated  and  the  approximate  results  are  within 
2%  error  of  the  exact  results  presented. 

The  parameters  ffhv/^hh  P  essentiaJly  unchanged  when  the  fractional 
volume  of  the  scatterers  changes  from  /  =  0.5%  (Table  1)  to  /  =  1%  (Table  2). 
Results  for  /  up  to  a  few  percent  have  also  been  evaluated  and  are  eJso  within 
a  few  percent  difference  with  the  numbers  shown  in  Table  1  for  those  two  pa- 
rauneters.  This  indicates  that  P  serve  as  useful  parameters  in 

identifying  and  distinguishing  the  geometry  of  scatterers  for  the  remote  sensing 
of  vegetation  where  /  is  usually  smaller  than  a  few  percent. 

APPENDIX  A;  EFFECTIVE  PERMITTIVITY  UNDER  THE  LOW  FREQUENCY 
APPROXIMATION 

Assuming 

|((e//)t;|  1  (Al) 

then,  from  (16), 

hff  -  +  (A2) 

Substituting  (A2)  and  (19)  into  (15),  then 

{D{f))  =  ,,7{E(f))  +  kl€o  (A3) 

Assuming  the  frequency  is  low  enough  so  that  the  variation  of  mean  field  (E[f')) 
is  negligible  within  a  few  correlation  lengths,  we  can  approximate  (A3)  as 

(Bif))  =  +  klto  J  (tir)  •  PSBg[fy)  ■  ?(r')))  •  {E(f))  (A4) 

Therefore,  we  can  define  the  effective  permittivity  of  the  medium  under  the  low 
frequency  approximation  as 

-  =  =(0) 

<«//  -  +  <0  (eff  (A5) 

where 

S//  dr'  ({(r)  .  (r,f')  •  f(f'))  (A6) 


APPENDIX  B:  EXPRESSION  FOR  and  5, 

For  the  correlation  function  (37)  in  the  local  coordinate,  the  integration  of  a  more 
general  case  than  (27)  has  been  carried  out  analytically  by  Nghiem  [8].  In  this 
paper,  results  are  presented  for  isotropic  background  medium  only  and  are  given 
by 


/s=  - 


+  (7’  -  IK?  -  7/^ 
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For  prolate  spheroids  (c  >  a) 

Ac 


1 


2e  +  In 


1  -  e 

l+e 


with 


and  Aq  can  be  obtained  from  the  following  formula, 

2 

2Aa  -}-  Ac  —  A 
a^c 
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(C4) 

(C75) 

(C6) 
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1.1  Introduction 

With  advances  in  polarimetric  radar  technology,  polarimetry  has 
become  inqKirtant  to  the  remote  sensing  of  geophysical  media.  IVilly 
P<datimetric  radar  signals  convey  additional  information  regarding  the 
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